Elastic deformation due to tangential capillary forces 
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A sessile liquid drop can deform the substrate on which it rests if the solid is sufficiently "soft" . fn 
this paper we compute the detailed spatial structure of the capillary forces exerted by the drop on 
the solid substrate using a model based on Density Functional Theory. We show that, in addition 
to the normal forces, the drop exerts a previously unaccounted tangential force. The resultant effect 
on the solid is a pulling force near the contact line directed towards the interior of the drop, i.e. not 
along the interface. The resulting elastic deformations of the solid are worked out and illustrate the 
importance of the tangential forces. 



I. INTRODUCTION 

A drop of liquid can adhere to a solid surface due to 
molecular interactions. The microscopic structure of the 
forces exerted on the liquid is accounted for by classical 
wetting theory, e.g., using the concept of disjoining pres- 
sure From a macroscopic perspective this induces a 
normal force localized near the contact line of magnitude 
7 sin 0, 0] J where 7 is the liquid- vapor surface tension 
and 9 the contact angle. This attractive force (per unit 
line) is compensated for by a repulsive Laplace pressure 
spread out over the contact area of the drop, giving a zero 
resultant force. Experimental evidence for such nontriv- 
ial spatial structure of capillary forces can be inferred 
from the elastic deformation of the solid below a liquid 
drop [1-[T^. Figure [T^ (taken from Q), shows that clas- 
tic deformation is small, but measurable, on sufficiently 
soft substrates. 

One should note, however, that elastic deformation is 
induced by forces exerted on the solid. This means that, a 
priori, one cannot apply the classical theories of wetting, 
which account for the forces exerted on the liquid. The 
goal of the present paper is therefore to extend the stan- 
dard formalism of wetting to compute the microscopic 
forces exerted on the solid. It will turn out that, in- 
deed, there are some crucial differences when changing 
the "system" to which the forces are applied from liquid 
to solid. We will work out the consequences for the re- 
sulting elastic deformations, which we show to be funda- 
mentally different from existing theoretical predictions. 

There has been an intrinsic interest in the elasto- 
capillary interactions such as drops on a soft surface 
[l7[. For example, the contact angle of the drop on a 
soft surface is slightly different from Young's equilibrium 
angle 6 [l8l - l2l| . The energy balance now has an elastic 
contribution and for a small deviation from Young's an- 
gle A9 one finds A9 ^ ^/{ER), where E is the Young 
elastic modulus of the system and R the radius of the 
drop. Other effects are enhancement of contact angle hys- 
teresis 12 21 . visco-elastic dissipation during drop spread- 
ing 1^ 2^ or enhanced nucleation density in condensa- 
tion processes [2^ . To estimate the significance of elastic 



deformations it is instructive to consider the characteris- 
tic length scale £ = 7/i?, from the ratio of surface tension 
and the elastic modulus, which determines the scale of de- 
formations. For polycrystalline materials, for which the 
elastic modulus is determined by the interaction between 
molecules [E ~ GPa), such elastic deformations are irrel- 
evant {I ~ 10""'^^ m). However, for solids whose elasticity 
is of entropic origin (rubbers or gels with E < IQ kPa) 
deformations extend into the micron range and have in- 
deed been observed experimentally [1, H, [ll|, [H, [13, HI] . 

Interestingly, the elastic problem cannot be solved 
within a macroscopic framework where the effect of sur- 
face tension is represented by a localized force. This is 
because the elastic response to a Dirac (5-function gives 
an infinite displacement at the contact line. This sin- 
gularity can be regularized by introducing microscopic 
physics, cither in the solid phase p7| or in the liquid 
phase. The simplest and most common approach is to 
spread out the pull of the contact line over a small re- 
gion of finite thickness w. By assuming that this force 
is uniformly distributed in this thin region, Rusanov Q 
and Lester Q have solved the elastic displacement pro- 
file analytically. The resulting scale of displacement near 
the contact line becomes ^ £ ln(i?/w), illustrating the 
singularity as w — > 0. As pointed out by White [J|, the 
precise distribution of capillary forces can be obtained 
explicitly when introducing a microscopic model for the 
molecular interactions. It was argued that the normal 
traction is simply equal to the disjoining pressure Ii{h). 
This indeed gives a pressure peak localized near the con- 
tact line, whose integral is equal to ^sm.9. Once more 
it should be emphasized, however, that this involves the 
forces exerted on the liquid, and not the forces exerted on 
the solid, as required to obtain the elastic deformation. 

In this paper we compute the forces exerted by the 
liquid on the solid using Density Functional Theory in 
the sharp-kink approximation |28j . This approach has 



the merit that (for small contact angles) all results can 
be expressed in terms of the disjoining pressure and can 
thus be related to the usual wetting theory. Our central 
finding is that, besides a normal component, there is also 
a nonzero tangential component exerted on the solid, as 
depicted in Fig.lT}). The physical origin of this tangential 
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FIG. 1. (a) Shape of an ionic liquid drop on a soft PDMS 
surface, determined experimentally by Pericet-C'amara et. al. 
01 using laser scanning microscopy. The aspect ratio has been 
exaggerated by a factor 15. (b) Macroscopic representation 
of the capillary forces below a small liquid drop and resulting 
elastic displacements. A localized normal force of magnitude 
ysinO pulls along the contact line and is balanced by the 
Laplace pressure, which is applied to the wetted surface area. 
In this paper we show that there is also a nonzero tangential 
force of magnitude 7(1 + cos6'). The resulting force points 
towards the interior of the liquid, at an angle a to the hori- 
zontal. The sum of these forces induces a small deformation 
of the interface, (c) Physical origin of the tangential force. 
Solid molecules just below the contact line experience a resul- 
tant attraction towards the interior of the drop, where there 
is more liquid. This effect is particularly striking when the 
contact angle 6 is above tt/2. 

force can be understood by considering solid molecules 
directly below the contact line: there is always more liq- 
uid present at the interior of the contact line than at the 
exterior, yielding an "inward" attractive interaction with 
a nonzero tangential component. Near the contact line, 
this leads to the following normal (F„) and tangential 
(Ft) force components: 

K/L = 7sin0, (1) 
Ft/i = 7(l + cos0), (2) 



where L is the length of the contact line. Hence, the re- 
sultant force points towards the liquid (Fig. [TJ;) and is 
not directed along the interface. The angle that char- 
acterizes the direction of the force follows from tan a = 
sin6'/(l -I- COS0), which simply gives a = 9/2. This force 
was not considered before in the literature and we there- 
fore investigate the way it affects the elastic deformation 
of the substrate. 

It is worth emphasizing that the capillary forces pre- 
dicted by the Density Functional Theory model are per- 
fectly consistent with thermodynamics and do not lead 
to a violation of Youngs law. One should bear in mind 
that (1,2) represent the forces exerted on the solid by 
contrast, Youngs law reflects the equilibrium shape of 
the dcformable liquid and thus involves forces exerted on 
the liquid 2^. In fact, the same model was previously 
shown to yield Youngs an gle w hen computing the equi- 
librium shape of the liquid (30l. l3l| and thus yields a con- 
sistent thermodynamic picture. Let us further note that 
the presented results are fundamentally different from 
the tangential contribution described in Refs. [1, 0, Q. 
Those refer to forces originating from the difference with 
Young's angle induced by deformation itself. This is a 
higher order effect with respect to the basic tangential 
capillary forces and are smaller by orders of magnitude 
(^ i/R). 

The paper is organized as follows. Section HIl develops 
the formalism used to calculate the forces on the solid, 
and points out the relation with disjoining pressure and 
thermodynamic pressure. In Sec. IIIII we then quantify 
the tangential and normal forces, leading to (|ll2p . and 
illustrate the results for a Van der Waals model for the 
interactions in Sec. IIVI In Sec. |V] we then show how the 
tangential capillary force affects the elastic deformation 
of the solid in the limit of small deformations, and the 
paper closes with a discussion. 

II. DENSITY FUNCTIONAL THEORY IN THE 
SHARP-KINK APPROXIMATION 

Here we introduce the formalism needed to compute 
the forces exerted by a liquid drop on a solid. This is 
based on Density Functional Theory in the sharp-kink 
approximation, in the same spirit as f28l . [30l - l33| . We will 
define four important pressures: the disjoining pressure 
n, the thermodynamic pressure Pl, the repulsion pr and 
the Laplace pressure. These will be used in Sec. IIIII to 
make explicit predictions for the capillary forces on the 
solid. 



A. Repulsion pr 

From standard continuum mechanics one can express 
the force F that the liquid exerts on the solid as 

F = - / drVf^is - / dAp„ (3) 
Js Jas 
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where S and dS are the volume and boundary of the 
sohd domain respectively. Long-range (attractive) inter- 
action between liquid and solid molecules gives rise to 
a potential (pLS inside the solid i.e. an energy per unit 
volume. The short-range (repulsive) interactions are as- 
sumed to lead to isotropic stresses and to determine com- 
pletely the pair correlation functions at a given density. 
The liquid therefore exerts a contact pressure Pr on the 
solid. Such an explicit separation in long-range attrac- 
tion and short-range repulsion is a standard approach 
in Density Functional Theory 3J], and forms the basis 
for our analysis. We furthermore use the so-called sharp 
interface approximation [28l . Isol . [s^ , where the different 
phases are assumed to be homogeneous and separated by 
an infinitely thin interface. In this framework, one can 
express the liquid-on-solid potential as 



"LSI 



(r) ^ PlPs / drVLs(lr-r'l), 



(4) 



where (pLS is the effective solid- liquid interaction poten- 
tial, which takes into account the pair correlation func- 
tion (throughout we assume the interaction to depend 
only on the intermolecular distance |r — r'|). and ps 
are the liquid and solid densities respectively. The inte- 
gral runs over the entire liquid domain £, which implies 
that (j>LS is a functional of C. More generally, one can 
define 



{y)^PaPi3 / drVa/3(|r - r'l), 



(5) 



as the potential in phase /3 due to phase a. Note that in 
general 4>af} ^ 4>i3a, unless the domains of a and /3 have 
an identical shape. 

The calculation of the force F (or to be precise the nor- 
mal component Fn), still requires the repulsive pressure 
Pr- This can be obtained from the equilibrium condition 
inside the liquid. 



V(p^ 



) = 



(6) 



so that the tota.1 potential Pr + 4>ll + 4>sl must be ho- 
mogeneous [2^, [sol . Here it is assumed that the incom- 
pressibility of the liquid results from the repulsive term 
Pr, which adapts its value according to the local values of 
the liquid-liquid and solid- liquid interactions 4>LL,'t>SL- 



B. Thermodynamic pressure Pl 

The existence of a quantity conserved throughout the 
liquid (the total potential) allows to generalize the con- 
cept of thermodynamic pressure Pl of the liquid phase. 
Let us consider the case where the liquid is in equilibrium 
with its vapour phase. The constant of integration in ^ 
can then be determined from the boundary conditions at 
the interfaces. In the present framework repulsion is de- 
scribed as a contact force, as opposed to attraction whose 
influence is spread overs a few molecular diameters. As 



a consequence the repulsive potential p^ must be contin- 
uous, while there is no such condition for (p^L and (j^SL- 
Inside the vapour the density is so small that interac- 
tions can be neglected compared to the kinetic pressure 
Py. Therefore, the only force transmission by the vapor 
on the liquid is through pr and we simply find pr ~ Pv 
in the vapor. When the liquid-vapor interface is flat, we 
thus find Pr ~ Pv = Pl at the both sides of the inter- 
face. If in addition the liquid film is macroscopically thick 
(meaning that (/)sl = at the interface), we thus obtain 
inside the liquid 



Pr + (t)LL + (f>SL =Pl+ nLL(O), 



(7) 



where we denote IIll{0) as the value of (pLL at the edge 
of an infinite half space of liquid. For a detailed definition 
of IIq,^ we refer to the subsequent paragraphs. 

We propose to take ([7]) as a definition of the thermo- 
dynamic pressure Pl in the liquid. It is, as expected, a 
characteristic of the phase and it can be identified to the 
true thermodynamic pressure in the macroscopic limit. 



C. Disjoining pressure 11 

The disjoining pressure is defined as the pressure acting 
on a flat liquid film of thickness h, due to the influence 
of the solid substrate (Fig. It can be related to the 
potential Ilai3{h) at a distance h from an infinite half 
space: 



n„0(/i) 



dx 



dy 



2tt 



( 1 



dz <Pa/3(|r|) 



(Pai3{r)dr . 



(8) 



Contrarily to 4>ai3 7^ 'Ppai recover the symmetry prop- 
erty Ila^ = n^Q,. Within the sharp-kink approximation 
discussed above, the disjoining pressure n(ft,) can then 
be computed as [sil. l32j 



n(/i) = ^5L(/l)-^LL(M• 



(9) 



From this definition it is clear that n(/i) quantifles the 
change in energy when replacing part of the liquid by 
solid molecules. In other words, it is the correction to 
the total potential at the liquid vapour interface due to 
the presence of the solid at a distance h. 

To relate to the macroscopic concepts of surface ten- 
sions and contact angles, one can use the following inte- 
gral properties [3l| : 



/•oo 

/ (/i) d/i = -27, 

^0 



HsL {h) dh 



-isv + ISL - 7 
-7 (1 + cosO) . 



(10) 



(11) 



Note that the solid-liquid interaction only involves the 
difference Jsl — isv, and not the surface tensions indi- 
vidually. Here we made use of Young's law, which was 



4 



vapour „ _ p 




vapour 



liquid 



oo 



Pr = Py+^fK-Iisd'^) 



X 




FIG. 2. Repulsion Pr for a flat liquid film of thickness h. The 
potentials Ull and Hsl are defined in the text. 




h 



h 



T 



FIG. 3. The potential energy at the top surface of a liquid 
film of thickness h can be decomposed into three contribu- 
tions. The solid-liquid interaction (j)sL is simply the potential 
due to an infinite half-space of solid located at distance h. 
The liquid-liquid interaction (fj^L can be expressed into two 
contributions of semi- infinite liquid phases: in one case it is 
located at a distance h, while in other case it is at a zero 
distance. 



explicitly validated in the context of the present DFT 
model (sol . [3ll |. Combining the two equations, one ob- 
tains the usual normalization of the disjoining pressure: 



n [h) dh 



-isv + ISL + 7 
7 (1 — cos 61) . 



(12) 



Note that, in the case of a description that includes a 
precursor film of thickness /iq, the lower limit of integra- 
tion should be ho instead of for each of the integrals. 

We can now provide the detailed structure of the po- 
tential fields inside the flat liquid film (Fig. [5]). Once 
again, the density inside the vapor can be neglected as 
compared to the kinetic pressure and we simply find 
Pr = Pv Crossing the liquid vapour interface, the repul- 
sive pressure is continuous. The attractive potential at 
the liquid- vapour interface can be decomposed as shown 
in Fig. [31 The attraction by the solid leads to a potential 
4>sL = {h). Namely, the solid represents an infi- 
nite half-space at a distance h. The self-attraction by 
the liquid film can be also expressed as a function of the 
interaction induced by half a space of liquid located at 
a distance h: from the geometric construction in Fig. [3] 
one finds = Ull (0) — Hll {h). We finally get the 



FIG. 4. Repulsion Pr for a macroscopic film of liquid that has 
a curved liquid-vapour interface. By convention, the curva- 
ture K is positive on the schematic. 



expression of the thermodynamic pressure from Eq. ([7]) 

Pl = p, + nsL(/i)-nLL(/i) 



(13) 



where we used the definition of the disjoining pressure 
(PI . For a thin fiat film, the thermodynamic pressure in- 
side the liquid is thus different from the vapour pressure. 
The difference is precisely the disjoining pressure. This 
justifies and explains the concept. 



D. Laplace pressure 

When the interfaces are curved, it can be shown by 
expansion of the domain of integration [2^, [s^l that there 
is a deficit of energy that is exactly equal to the Laplace 
pressure. Consider the liquid film shown in Fig. 31 which 
is assumed to be much thicker than the molecular size. 
The liquid-vapor interface presents a curvature n. If the 
curvature is macroscopic, /io|k| 1, the potential due to 
the liquid film turns out to be 



')LL =nLL(0)+7K 



(14) 



If K > 0, the liquid domain is convex so that the volume of 
liquid attracting a given element of liquid at the surface is 
smaller than in the flat case. Therefore, \(j>LL\ < |nLL(0)| 
- note that (j)LL and 11^^ are negative. We finally get the 
expression of the thermodynamic pressure from Eq. ([71): 



(15) 



Indeed, one recognises the standard Laplace pressure. 
We see here that, just like the disjoining pressure, the 
Laplace pressure results from the geometry of the attrac- 
tive volume. 



E. Local approximation 

In principle the above equations are sufficient to evalu- 
ate the liquid-on-solid force F, once the molecular inter- 
action ipLs is specified. However, the technical difficulty 



is that both (p^s and Pr are functionals of the hquid do- 
main C. In turn, the shape of this domain has to be 
found self-consistently such that the hquid is in equihb- 
rium (28l . Isol - fssj i.e. by minimizing the total free energy. 
For a hquid-vapor interface, the equihbrium condition 
is that the interface is an iso-potentiaL As pr = Pv is 
constant, this imposes that 4>ll + ipSL is also constant 
along the free surface. This can readily be inferred from 
Eq. ([7]). It was shown in [13, [U that this framework 
indeed leads to Youngs law for the macroscopic contact 
angle of the liquid at equilibrium. 

Th equilibrium condition can be simplified for inter- 
faces of macroscopic curvatures and small slopes [3l|, [s^ . 
Then, the nonlocal functionals reduce to Junctions of the 
local thickness h and local curvature k. The two geomet- 
rical effects previously described add to each other: the 
capillary pressure is the sum of the Laplace pressure and 
the disjoining pressure. At the free surface, the solid 
induced potential reads (jisL = ^SL{h) and the liquid 
induced potential reads: 



^LL=^LL{0)-IiLL{h)+lK 



(16) 



In fact, one can show from a geometric construction that 
this local approximation is valid when h 3> h^h' ^ which 
for small slopes spans the entire liquid domain. 

One finally obtains the thermodynamic pressure Pl 
from Eq. (O: 

Pl=Pv- HllW + ^SL{h) + in 

= p„ + 7K + n(/i) . (17) 

This equation can be used to compute the equilibrium 
shape of the hquid domain, expressed as h(x,y). This 
will be done explicitly in Sec. IIVI where we compute the 
shape of drops and capillary forces for a specific choice 
ofU{h). 



III. CAPILLARY FORCES ON A SOLID 

We will now explicitly compute the force exerted by 
the liquid on the solid, F, starting from The Density 
Functional model discussed above will allow for direct 
evaluation of F, and also provides the detailed spatial 
distribution of the capillary forces in the vicinity of a 
contact line. Throughout this section we assume that 
the solid is undeformed. 
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FIG. 5. Interaction of a volume element within the solid (at 
position 2 < 0) with the liquid film of thickness h. The energy 
(i>LS can be decomposed into two potentials due to a semi- 
infinite liquid phase: in one case the solid element is at a 
distance —z from the semi-infinite liquid phase, whereas in 
the other case it is at a distance h — z. 
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FIG. 6. (a) For a fiat film the force exerted by the liquid 
on a volume element of solid is in the direction normal to 
the interface, (b) When the liquid-vapour interface is not 
parallel to the substrate, the left-right symmetry is broken a 
tangential component appears. The bias is towards the side 
where the amount of attracting liquid is greater. For a postive 
slope dh/dx, this is oriented in the positive x direction. 



4>LS then arises in the direction normal to the interface. 
Therefore, Eq. ([3|) can be integrated along the normal 
direction yielding a normal force inside the solid, 



dA [4)LS + Pr] 



(18) 



as 



where the integrand of ([T5)) has to be evaluated at the 
solid-liquid interface. Since the volume of liquid C is 
almost a semi-infinite domain bounded by a flat inter- 
face, we can write (ji^s — ^ls{^) — 115^(0) and (j)LL — 
nLL(O). Using Eq. 0, we obtain p^ = - 05^(0). 
Finally, this yields the total normal stress acting on the 
solid (Fig. g]): 



In 



-Pr 
-P, 



-Pl 



(19) 



which is the force per unit area. For a convex liquid inter- 
face (k > 0), the solid is submitted to a higher (Laplace) 
pressure. 



A. Capillary forces under a macroscopic film 

Far from a contact line, when the liquid film is macro- 
scopic, one expects the force per unit area acting on the 
solid to be simply the thermodynamic pressure Pl , trans- 
mitted at the interface. Here we check that this property 
is verified by the Density Functional Theory model. 

We consider again a perfectly flat solid-liquid interface 
separating an infinite half space of solid from a macro- 
scopic drop of liquid (Fig. 0]). The only gradient of 



B. Tangential capillary force near a contact line 

To identify the capillary forces near a contact line we 
now consider a microscopic liquid film. We first derive 
the component of the liquid-on-solid force that is tangen- 
tial to a flat solid interface. From ^ one can infer that it 
results from the long-range term involving ^ls, since tr- 
acts along the direction normal to the surface. Suppose 
that the contact line extends in the y-direction, so that 
the liquid- vapor interface can be written as h(x). Since 
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d/dy = 0, the tangential force Ft is along the a;-direction 
and can be expressed as 



Ft 



dy 



dx I dz— [(j>Ls] , (20) 



where the integrals run over the entire solid domain. 

We now apply the local approximation to determine 
(/)LS- At a position z < inside the solid, the potential 
due to the finite slab of thickness can be expressed as 
(Fig. EI: 



t>LS = IIsL (-z) - UsL {h {x) - z) 



Here we used the symmetry property IIq,^ 
sequently, 



n 



d_ 

dx 



(21) 



Con- 



(22) 



where h' = dh/dx. This implies that there will be a finite 
contribution of the tangential force only when the liquid 
film thickness varies with x {h' ^ 0). This breaking of 
left-right symmetry, illustrated in Fig. [51 induces a tan- 
gential force towards the side where the greater amount 
of liquid is present. Replacing ((22|) in ((20)) . we can inte- 
grate out the 2-direction and write the tangential force 



as. 



Ft = - dy 



dxh'UsL (h), 



(23) 



where we used the property 115^(00) = 0. 

From this expression one can draw two important con- 
clusions. First, one can interpret the argument of the 
integral as the tangential traction ft, defined as the tan- 
gential force per unit area: 



ft{x) = -h'nsL {h{x)). 



(24) 



This tangential capillary traction will be used below as 
input to compute the elastic deformation. Strictly speak- 
ing, this force is not localized at the solid-liquid surface, 
but it is spread out over a thin region in the 2-direction. 
In the clastic calculation we will treat ft as localized on 
the surface. Second, the total tangential force can be 
quantified exactly from Writing J dx h' — J dh, 

one can use the integral property (ITT|) . to find the cen- 
tral result Ft/L = 7(1 -I- cos 6*) (with J dy = L), already 
presented in Eq. 

Actually, Eq. ^ does not rely on the local approxi- 
mation and can be derived exactly, for arbitrary values 
of the contact angle 9. Starting again from Eq. ((20)) . the 
integration over x leads to: 



Ft = 



dy 



dz I 



(25) 



We assume that the liquid bulk is located in the pos- 
itive X direction, giving (t>Ls\x=oo = ^ls{~z)- At the 
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FIG. 7. The interaction of a volume element of liquid, located 
at the solid-liquid interface, can be decomposed into the inter- 
action exerted by the solid (/>s l , plus two interactions between 
the same liquid element with a semi-infinite liquid phase (in 
one case at a distance h and the other case at zero distance). 



same time, the dry solid or the precursor film is located 
in the negative x direction, h{—oo) = ho. Using the 
construction of Fig. [5] with film thickness ho, this gives 
4>Ls\x=-oo = TIls{-z) - TlLs{ho 
expressions in (PS)) . we get 



z). Inserting these 



Ft/L 







dz UsL {ho - z) 

dzUsL (5) =7(1 



COS 61), (26) 
where we once more used the normalization (jllll . 

C. Normal capillary force near a contact line 

The normal force near a contact line can again be 
computed from (jl8[) . but now we take into account that 
the liquid thickness is finite and lies within the range 
of molecular interaction. For a liquid profile h{x), the 
expression for the normal force becomes 



Fn = - dy 



dx [pr + (t)Ls)^=o ■ 



(27) 



From Fig. [5l we compute the potential inside the solid, 
located at the solid-liquid interface as 

0LS = nis(O) - IlLs{h) = nsL(O) - IisL{h) . (28) 

To determine Pr we again use the equilibrium condition 
([7|) inside the liquid, and evaluate it at the solid-liquid 
interface. This requires the potentials (pLL and (psL at the 
solid- liquid interface, as sketched in Fig.[71 The potential 
due to the solid is simply 4>sl = ^sl{^), since the solid 
is an infinite half space. Using the local approximation, 
the liquid-liquid interaction can be expressed in terms of 
Ull (cf. Fig.©: 



<^LL-nii(o)-nii(/i) . 

Hence, one obtains from ([7]): 

= PL+nLi(/l)-HsL(0). 



(29) 
(30) 



7 



This gives the final expression for the normal force: 
Fn^- fdy f dx{PL + TlLLih) - lisUh)} 

= -fdyf dx{PL~nih)}. (31) 

J J —OC 

This result can be interpreted as follows. The argu- 
ment of the integral represents the normal traction /„ 
that can be used to compute the elastic deformation: 



fnix)=-PL+Tlih{x))^ 



(32) 



This is indeed identical to the normal traction proposed 
by White (3]. The disjoining pressure contribution n(/i) 
acts only in the vicinity of a contact line, where h is in 
the range of the molecular interactions. It is responsible 
for the "pulling" force at the contact line. The second 
term, P^, is simply the thermodynamic pressure in the 
liquid, which was already identified in Sec. IIII Al for a 
macroscopic drop. 

To quantify the magnitude of the pulling force, we now 
consider a two-dimensional profile for an infinite drop 
with zero external pressure, i.e. Pl = Pv — 0- Some nu- 
merical examples will be presented in the next section (cf. 
Fig. [5]). The equilibrium profile is then given by a bal- 
ance between Laplace pressure and disjoining pressure, 
which according to (|17p can be written as 



n{h) + "fK n(ft.) -7 

This can be integrated to Q 
h' 

7 



h" 



(1 + /l'2)-^/2 



= 0. 



(33) 



(i + /,'2)i/2y^^ 



dxU{h), (34) 



where we used h'(—oo) ~ corresponding to the dry 
solid or the precursor film. Thus with h'{oo) = tan^, 
the left-hand side becomes 7 sin 6*. The right-hand side 
is an integral over the normal traction /„. Hence, the 
integrated normal force at the contact line becomes, for 
Pl = 0, 



Fn ^ J dy J dxll{h{x)) = J dy ^sinl 



(35) 



which is the anticipated result ((T]). 



IV. NUMERICAL RESULTS 

A. Drop shape and capillary traction 

From the analysis above it is evident that in order to 
compute the spatial structure of tangential and normal 
tractions, one needs to know the exact shape of the drop. 
This can be done once the disjoining pressure has been 
specified. As an illustration of our results we consider a 



disjoining pressure n(/t) that corresponds to a Lcnnard- 
Jones (9-3) potential fsa. [36j. 



TlsL (h) 



A 
ho 

TIll (h) = ^ 
ho 



Il{h) 



B-A 
ho 



(36) 
(37) 
(38) 



without any correction due to the pair correlation func- 
tion. The molecular length ho is the equilibrium thick- 
ness of the precursor film and sets the length scale for 
the interactions. For simplicity we take the same func- 
tional form for liquid-liquid and solid-liquid interactions, 
but this is not required. The second term is the famil- 
iar long-range van der Waals attraction, while the first 
represents the short-range repulsion. Such a repulsion 
is needed to prevent matter from collapsing and ensures 
that the surface tensions remain finite, cf. (|10lllll2p . In 
this case, the repulsive term balances the attraction at a 
molecular length ho, which selects the equilibrium thick- 
ness of a precursor film. To illustrate the robustness of 
our results. Appendix El also discusses a disjoining pres- 
sure for which the repulsion is modeled by a vanishing 
pair correlation function. Despite the fact this model 
does not lead to the formation of a precursor film, the 
resulting tractions on the solid turn out very similar. 

The constants A and B are determined by the integral 
properties (|10lllll2p . as A = |7(l-hcos6') and B = 
-y7. Below we compute the shape of the drop and the 
corresponding capillary traction /„.t for two cases: an 
infinite two-dimensional drop and a finite axisymmctric 
drop. 



B. Infinite two-dimensional drop 

The drop shape for an infinite two-dimensional drop 
can be determined by the equilibrium condition (j33p . 
The boundary condition used for the solution of the drop 
profile is that the interface asymptotically joins the pre- 
cursor film, /i(— 00) — ho- This automatically gives the 
correct contact angle as h'{oo) = tanO. 

Figure [8] shows the drop shape and the capillary trac- 
tion for two contact angles, 9 = 20° and 70°. The solid 
line (green) represents the shape of the liquid h(x) nor- 
malized by the precursor film thickness ho- The peaked 
curves are the traction /„ (dash-dotted line in blue) and 
ft (dashed line in red), normalized by the natural pres- 
sure scale in the model fo = J /ho- Indeed, one can see 
that the capillary forces are localized in the vicinity of the 
contact line, in between the precursor film and the macro- 
scopic drop. The width of the peaks is set by the length 
scale ho- The most important result highlighted in these 
figures is that for a small contact angle, the tangential 
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traction is significantly larger than the normal traction - 
this is consistent with Eqs. (|ll2p since the total force (per 
unit length) is obtained by integrating over the peaks. In 
fact, this integral property is determined by the values of 
the surface tensions and holds irrespective of the choice of 
disjoining pressure. Only the detailed shape of the peaks 
is affected by the choice of disjoining pressure - compare 
e.g. to the example given in Appendix]^ It is worth not- 
ing that, although the local approximation may not be 
fully quantitative, the figure is still representative since 
the localization of the capillary forces and their relative 
strengths are independent of this approximation. 
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FIG. 8. (Color Online) Drop shape and capillary tractions for 
an infinite drop with (a) 9 = 20° and (b) = 70°. In both 
plots the drop profile h{x)/ho is represented by a solid line 
(green), the tangential traction ft/ fo by a dashed line (red) 
and the normal traction fn/fo by a dash-dotted line (blue). 
Tractions are normalized by the pressure scale fo = j/ho- 



C. Axisymmetric drop 

For three-dimensional axisymmetric drops, the equilib- 
rium equation (|17p for the drop shape h{r) reads [37| 

jK + U{h) = PL-Pv, (39) 
where the expression of the curvature k is now 
h" h' 



(1 + /i'2)^/2 



(40) 



This equation is solved numerically with the boundary 
conditions (^')r=o ~ ^ (symmetry) and (/i')r->-oo ~^ 
(matching onto a precursor film). The same Lennard- 
Jones disjoining pressure is used. Different drop volumes 
can be obtained by varying — Py 



Figure shows the drop shape h{r) and normalized 
tractions ft and /„ for a contact angle 9 = 20°. Simi- 
lar results are obtained for different volumes and contact 
angles. The thickness of the precursor film is slightly dif- 
ferent from ho due to the finite volume, and is set by 
n (h) ~ Pl — Pv ■ Once more the tractions arc strongly 
peaked near the contact line. The negative sign for /( 
signals that it acts in the negative radial direction, i.e. 
towards the centre of the drop. A key difference with 
respect to infinite drop of Fig. [S] is that the normal trac- 
tion attains a nonzero constant value at the interior of 
the drop. This is simply the Laplace pressure (equal to 
0.00617/ft.o here), with a minus sign. 

This figure forms the quantitative basis for "cartoon" 
of the capillary forces that was presented in Fig. [T] From 
the sign of the normal traction, it is clear that inside 
the drop the solid is pushed vertically downwards by the 
Laplace pressure, whereas at the contact line it is pulled 
upwards by a force that decays rapidly on either sides of 
the contact line. Similarly, there is a localized tangential 
force pointing towards the interior of the drop. 



V. ELASTIC DEFORMATIONS BELOW A 
LIQUID DROP 

With the drop shapes and the tractions on the solid 
in hand, one can compute the elastic response of the 
solid substrate. We will use the tractions ft and /„ from 
Eqs. ([M)) and as obtained for the undeformed solid. 
Hence, we assume that the deformation of the solid does 
not affect the resulting capillary traction, which in prac- 
tice means j/ho ^ E. In addition, we treat the tractions 
as perfectly localized at the interface. Ideally, one would 
like to treat the solid as a semi-infinite elastic body and 
consider the tractions as a simple load per unit length. 
This cannot be done, however, since such a line-loading 
gives rise to logarithmically diverging displacements in 
the far field [38| . This problem can be avoided by con- 
sidering a finite amount of liquid, which is commonly 
done using axisymmetric drops of finite volume (39| . An 
alternative is to consider an elastic medium that has fi- 
nite thickness Q, which is mathematically much more 
involved. As already mentioned, none of this literature 
considers the elastic deformations due to tangential cap- 
illary traction. 

Here we consider elastic deformations induced by ax- 
isymmetric drops, whose tractions are functions of r only. 
The elastic response then follows from standard expres- 
sions for semi-infinite elastic bodies that can be found 
e.g. in [3^. We distinguish the following cases: 

• normal displacement due to normal traction (Sn.n) 

• normal displacement due to tangential traction 

• tangential displacement due to normal traction 



9 




50 r/ho 100 



4 

-5 

5t,t 
-10 







(c) 











50 



r//io 



100 



150 



FIG. 9. (Color Online) Radial profiles of height, stresses 
and displacements for an axisymmetric drop, (a) Drop profle 
(solid green line); tangential traction ft (dashed red line); 
normal traction /„ (dashed-dotted blue line), (b) Normal de- 
formation due to normal (dashed-dotted blue line) and tan- 
gential (dashed red line) tractions, (c) Tangential deforma- 
tions due to normal (dashed-dotted blue line) and tangential 
(dashed red line) tractions. 



• tangential displacement due to tangential traction 
(^m) 

The normal displacement at the solid surface due to 
normal traction can be obtained as follows: 

{r) = —^l —^K{m)f^{r')dr' (41) 



■kE Jq r + r' 

Here v is the Poisson ratio, E is the Young's modulus 
of the material and K ijn) is the elliptic integral of the 
first kind with m = j^^^^py- Similarly, the tangential 

displacement due to tangential traction can be computed 

as 

4r' 



1 



ttE 
2 
m 



r + r' 

1 ) K{m) 



2 



Y{m) 



dr'. (42) 



where Y{rn) is the elliptic integral of the second kind. It 
should be noted that the Kernels display a logarithmic 



singularity at r = r', i.e. m = 1. This means that ap- 
proximating the tractions as Dirac 5 functions perfectly 
localized at the contact line would give a diverging de- 
formation at the contact line. Instead, with the tractions 
spread out over a finite width of order /iq, the maximum 
displacement scales as '-^ ln(i?//io)j where R is the drop 
size. 

By contrast the normal displacement due to tangential 
forces, 5n,ti can be obtained in a much simpler form, 
without invoking the elliptic integrals: 



5n,t (r) 



ttE 



ft{r')dr', (43) 



and similarly, the deformation St.n follows as 
(1 - 2;^) (1 + I') 



5t,n (r) 



E 



-fn{r')dr'. (44) 



These expressions no longer display a singular Kernel. 
Note that the prefactors of these two expressions vanish 
for V = 1/2, which corresponds to incompressible elastic 
materials. 

From these expressions one can identify a typical 
length scale le for the clastic deformations. Since in- 
tegrals over /„ and ft arc proportional to the surface 
tension 7, we define the elastic length scale as: 



(1-^^)7 



ttE 



(45) 



Based on this length scale, we can now tabulate the scal- 
ing of the deformations in different directions (Table |l| . 
In these expressions we actually account for the 9 de- 
pendencies of the normal/tangential tractions, according 
to Eqs. (HI) and This clearly shows the importance 
of the tangential capillary forces. For normal displace- 
ments, the effect of the tangential traction is comparable 
to that of the normal traction, in particular when the 
contact angle is small. For tangential displacements, 
however, tangential forces even dominate the effect of 
normal forces and omitting them induces a significant 
quantitative error. 



TABLE I. Scaling of the displacements in normal and tan- 
gential directions due to normal and tangential tractions, ex- 

(l-^'')7 



pressed in the elastic length = 









Normal traction 


sin 6 


i^TTsinS 


Tangential traction 


1^{1 + COS0) 


1-f cos 6* 



To further illustrate this, we have computed the de- 
formation profiles for 9 = 20°, v ~ 0.35, and for the 
disjoining pressure given by Eq. (j38p . For the numeri- 
cal evaluations of the elliptic integrals (necessary for ^„ „ 
and 5t^t)^ we used a method discussed in Johnson |38j 
that explicitly avoids the singular Kernel. The resulting 
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displacement profiles are shown in Figs.|9lDC. Normal dis- 
placement produced by tangential traction is comparable 
to that produced by the normal traction, whereas the 
tangential deformation produced by the tangential trac- 
tion is significantly larger than the corresponding contri- 
bution of the normal traction. The displacements Sn^n 
and 6t,t display a maximum near the contact line, whose 
magnitude depends on the drop size as ^ ln(i?/ft,o)- For 
micron-sized or larger drops, the maximum will thus be 
much more pronounced than shown here since Hq is of 
molecular scale. A peculiarity of S^^t and Sf ^ is that 
the displacements are nonzero only at the interior of the 
drop, and can thus not be measured outside the drop. 



VI. SUMMARY AND PERSPECTIVE 

In this paper, we have used a nanoscopic yet contin- 
uum approach - the Density Functional Theory in the 
sharp kink approximation - to derive the capillary forces 
exerted on a solid substrate by a liquid drop resting on 
it. The drop exerts a large tangential capillary force in 
the vicinity of the contact line, with a resultant magni- 
tude of 7 (1 -I- COS0) per unit length. Such a force has not 
been reported before and its physical origin is sketched 
in Fig. [TJ;: there is a resultant attraction towards the in- 
terior of the drop, where most of the liquid is situated. 
Clearly, this bias is of geometric origin and therefore does 
not rely on details of the sharp kink approximation. We 
thus expect that a more complete theory only gives minor 
quantitative differences. Let us note that our model re- 
covers the normal component of the capillary force, which 
is known to have a resultant jsinO. In the bulk of the 
drop, the Laplace pressure exerts a vertical push on the 
substrate. 

This novel picture of capillary forces is next employed 
to determine the elastic deformations of the solid. This 
is relevant when the substrate is sufficiently soft - i.e. 
when the elastic modulus is sufficiently small - to lead to 
nanoscopic or microscopic displacements. The principal 
limitation of our elastic calculation is that it only applies 
for small deformations, i.e., 7//10 ^ E. At larger de- 
formation, the resulting change in geometry of the solid 
will have a significant effect on the capillary tractions, 
as demonstrated in the recent work by Jerison et al [13 ■ 
Another approximation is that we treated the tractions 
as localized at the surface of the substrate. However, 
the corrections induced by the finite range of the inter- 
actions will have an effect only in a small region around 
the contact line, where e.g. the maximum displacements 
arc known to depend logarithmically on details of the 
regularization 0, I^, [s^ . 

The predicted tangential force is substantial and will 
have a measurable effect on the elastic deformations. 
Typically, experiments studying the solid deformation 
below a liquid drop measure the normal displacements 
[13, [m - see e.g. Fig.[T^. This is equivalent to S 
in the present analysis. These substrates often have a 



Poisson ratio close to u = 0.5, in which case the con- 
tribution due to the tangential force Sn,t is expected to 
be relatively small. By contrast, the tangential deforma- 
tions of the free surface are very sensitive to the tangen- 
tial force, as can be inferred from Fig. |9] and Table ID 
We thus suggest that, in principle, tangential displace- 
ments should provide a good experimental tool to access 
tangential capillary forces. 

Finally, let us emphasize that a tangential force be- 
low the contact line of magnitude 7 (1 -I- cos6') provides 
a striking perspective on capillary forces. We first note 
that this does not lead to a violation of Youngs law for 
the contact angle: the computed tractions are exerted on 
the solid, while Youngs law reflects the forces exerted on 
the liquid. However, when projecting the surface tension 
along the liquid-vapor interface, the tangential compo- 
nent would be J cos 6 [iO] and not 7(l-f cos^). Simi- 
larly, when considering a plate plunged in a liquid bath, 
the virtual work principle dictates that the resultant force 
on the solid in fact is 7 cos 9 ^ . All these thermodynamic 
results are correct and do not contradict our present anal- 
ysis. Namely, 7 (1-1- cos 0) is the force per unit length 
exerted below the contact line - there are other forces ex- 
erted on the substrate at other places. For instance, the 
tangential forces around a drop balance each other. The 
case of the plunging plate is more subtle. We show in the 
Appendix [BF that a normal force appears on the solid 
in places where the substrate is curved. For the plunging 
plate, there are thus very strong forces acting on the bot- 
tom edges of the plate. Indeed, one can show that these 
ensure that the overall force on the solid is in accordance 
with the virtual work principle. Once again, we thus ex- 
pect direct experimental evidence of the tangential force 
7(1-1- cos 9), by measuring the tangential compression of 
a soft elastic plate in a liquid. 

Acknowledgements: We gratefully acknowledge dis- 
cussions with Joost Weijs and Kees Vernier. 

Appendix A: Tractions for a different interaction 
potential 

In this appendix we verify that the tractions shown in 
Fig.[H]are representative in the sense that results are only 
mildly affected by the choice of disjoining pressure. In 
particular, we will use an interaction model that does not 
lead to a precursor film, but has a well-defined contact 
line at which h = 0. For this, we consider an interaction 

representing a van der Waals attraction that is regular- 
ized at molecular scale ho. The regularization stems from 
the vanishing pair-correlation function within a distance 
/iQ. For numerical convenience we follow [2^ [s^l and 
use a simplified form that is easily integrated. From this 
potential we can derive the disjoining pressure from 
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r h \ 




- - tan"^ 






2 




1 



h/ha 



Employing (|X2|) in ^ and (HH, we get 



CSL 



(A2) 

(A3) 
(A4) 



and hence the total disjoining pressure becomes 



7(1 



ho 







- - tan~^ 




2 





h/ha 



1 



(V/io)2 
.(A5) 

We consider the shape of an infinite two-dimensional 
drop corresponding to this disjoining pressure. For this 
we numerically evaluate (|33p with the equilibrium bound- 
ary conditions at the contact line, namely ft.(0) = and 
/i'(0) = 0. Figure [TUb shows the resulting drop shape 
as well as the tractions from (|24p and ([5^ . These re- 
sults should be compared to Fig. [8^1. First, one notes 
that the new model does not lead to the formation of a 
precursor film. Yet, the resulting tractions display very 
similar characteristics as those obtained from the model 
with precursor film. The total force (per unit length) 
is obtained by integrating over the peaks and these are 
once more equal to 7sin0 and 7(1 + cos6'), for the nor- 
mal and the tangential components. Only the detailed 
shape of the peaks is affected by the choice of disjoining 
pressure, though for each of the cases the peaks have a 
typical width ~ /ig. 

A peculiar feature of the model without the precur- 
sor film is that the presence of the contact line seems to 
induce a discontinuity in the tractions. Namely, outside 
the drop, a; < 0, the tractions are strictly zero. Inside the 
drop, the normal traction is finite since lim^_^o+ n(ft.) = 
n(0) 7^ 0, and hence leads to a discontinuity. We show, 
however, that this discontinuity is an artifact of the lo- 
cal approximation of the DFT model. This approxima- 
tion implicitly assumes small spatial derivatives "d/dx" , 
a condition that is clearly violated near the contact line 
in case there is no precursor film. Using the methods de- 
scribed in [3^ , we were able to numerically solve the full 
DFT model for the interaction (|A1[) . By this we mean 
that we solve the drop shape and tractions from the func- 
tionals described in Sec. without invoking the local 
approximation. The results are presented in Fig. llOb . 
Indeed, the tractions are very similar to those predicted 
by the local approximation, the main difference being a 
smoothening of the discontinuity near a; = 0. 



Appendix B: Curvature induced force 

In this appendix, we compute the force exerted by the 
liquid on a solid of curvature depicted in Fig. 1111 From 
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FIG. 10. Drop shape and capillary tractions for a two- 
dimensional infinite drop for a model that does not lead to a 
precursor film. Molecular interactions given by (|A1|) . contact 
angle Q — 20° . (a) Numerical results obtained from the "local 
approximation" . (b) Numerical results obtained from the full 
DFT model, without invoking the local approximation. See 
text for details. Overall, the results are very similar to those 
of Fig. [8^. Note that the local approximation induces a dis- 
continuity in the normal traction, which is not present in the 
full model. 



the point of view of the liquid, the same interface has 
thus a curvature ~k. In this case the expression of the 
normal stress on the solid is slightly modified with respect 
to the flat surface, cf Eq. p^ . as it picks up a curvature 
contribution: 



fn = -Pr 



lisL (h) dh 



^ ~Pr - (j)LS - li'^ + C0s9)k . (Bl) 

Considering the sub-system shown in dark gray in Fig. 1111 
the curvature correction term comes from the reduced 
volume that is attracted to the liquid, as compared to the 
flat case. If the solid is convex, i.e. k > 0, this reduction 
of attraction effectively yields an increased pressure. 

As previously, repulsion can be expressed using the 
total potential. Using the analogue of equation ([T4|) for 
the solid-liquid interaction, one express the potentials: 

•pLL = 11^^(0) - 7 K 
1 



(t)sL = nsL(o) -I- -7(1 + cose*) K , 

<pLS = nLs(O) -^7(1 + COS 9) K . 

By symmetry, we also have 11^5(0) = 11^5(0). 
Eq. [71 we get the expression of the repulsion: 
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FIG. 11. Repulsion pr at the surface of a curved solid. By 
convention, the curvature k is positive on the schematic. 
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